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ABSTRACT 

We present SHARC-II 350^ni imaging of twelve 24^m-bright (F24^m > 0.8 mJy) Dust-Obscured 
Galaxies (DOGs) and CARMA 1mm imaging of a subset of 2 DOGs. These objects are selected 
from the Bootes field of the NO AO Deep Wide-Field Survey. Detections of 4 DOGs at 350/im imply 
infrared (IR) luminosities which are consistent to within a factor of 2 of expectations based on a 
warm dust spectral energy distribution (SED) scaled to the observed 24/im flux density. The 350/im 
upper limits for the 8 non-detected DOGs are consistent with both Mrk 231 and M82 (warm dust 
SEDs), but exclude cold dust (Arp 220) SEDs. The two DOGs targeted at 1mm were not detected 
in our CARMA observations, placing strong constraints on the dust temperature: T^ust > 35 — 60 K. 
Assuming these dust properties apply to the entire sample, we find dust masses of w 3 x 10^ Mq. In 
comparison to other dusty z ~ 2 galaxy populations such as sub- millimeter galaxies (SMGs) and other 
5pi!teer-selected high-redshift sources, this sample of DOGs has higher IR luminosities (2 x 10^'^ Lq vs. 
6 X 10^^ Lq for the other galaxy populations) that are driven by warmer dust temperatures (>35-60 K 
vs. ~30 K) and lower inferred dust masses (3 x 10^ Mq vs. 3 x 10^ Mq). Wide-field Herschel and 
SCUBA-2 surveys should be able to detect hundreds of these power-law dominated DOGs. We use 
existing Hubble Space Telescope and Spitzer/IKAC data to estimate stellar masses of these sources 
and find that the stellar to gas mass ratio may be higher in our 24/im-bright sample of DOGs than in 
SMGs and other 5piteer-selected sources. Although much larger sample sizes are needed to provide 
a definitive conclusion, the data are consistent with an evolutionary trend in which the formation of 
massive galaxies at z ~ 2 involves a sub-millimeter bright, cold-dust and star- formation dominated 
phase followed by a 24/im-bright, warm-dust and AGN-dominatcd phase. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — galaxies: high-redshift 
— submillimcter 



1. INTRODUCTION 

In the local universe, the most bolometrically lumi- 
nous galaxies are dominated by thermal emission from 
dust which absorbs ultra-violet (UV) and optical light 
and re-radiates it in the infrared (IR) (jSoifer et al.l 
I1986I) . While rare locally, these ultra- luminous IR 
galaxies (ULIRGs) are more common at high red- 
shift fe.g..lFrancesch ini ct al.l l200j] : lLe Floc'h et al.ll2005l : 
iPerez-Gonza lez etal. 2005). Studies combining the im- 
proved sensitivity in the IR of the Spitzer Space Telescope 
with wide-field ground-based optical imaging have iden- 
tified a subset of this z 2 ULIRG p opulation that 
is IR-bright but also optic ally fain t (Yan et al. '2004'; 

w 

[2008; Dcv et al 
and Fiore et al. 



Houck et al. 2005; Wccd man et al.l 2006b; Fiore et al.i 



In particular, Dev et al-l (|2008D 
present a simple and economi- 
cal method for selecting these systems based on i?-band 
and 24/im Multiband Imaging Photometer for Spitzer 
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(MIPS; IRieke et al.l[200l data. IDev et all (gOOp) em- 
ploy a color cut of i? — [24] > 14 (Vega magnitudes; 
KiF^{24^im)/F^{R) > 1000) to identify objects they 
call Dust-Obscured Galaxies (DOGs). Applied to the 
8.6 deg^ Bootes field of the NO AO Deep Wide-Field Sur- 
vey (NDWFS) that has uniform MIPS 24/im coverage for 
(24/im) > 0.3 mJy , this selection yields a sample of 
«2600 DOGs, or «302 deg-^. 

The extreme red colors and number density of the 
DOGs imply that they are undergoing a very luminous, 
short-lived phase of activity characterized by vigorous 
stellar bulge and nuclear black hole growth. Spectro- 
scopic redshifts determined for a sub-sample of DOGs 
using the Deep Imaging M ulti-Object Spectrograph 
(DEIMOS; iFaber et al.l I2003D and the Low Re solution 
Imaging Spectrometer (LRIS; lOke et aD 119951 ) on the 
telescopes of the W. M. Keck Observatory (59 DOGs), 
as we ll as the Infrared Spectrometer (IRS; IHouck et al.l 
|2004[) on Spitzer (47 DOGs) have revealed a redshift dis- 
tribution centere d on z « 2 with a dispersion of cr^ « 0.5 
(IDev et al.ll2008l) . 

While DOGs are rare, they are sufficiently luminous 
(« 90% of DOGs with spectroscopic redshifts have 
Liji > 10^^ Lq) that they may contribute up to one- 
quarter of the total IR luminosity density from all 
z 2 galaxies (and over half that from all ULIRGs 
at z ~ 2) and may be the progenitor s of the most lu- 
minous (~4L*) pre sent-day galaxies (|Dev et al.l 120081 : 
iBrodwin et al.l 120081 ). Thus far, the efforts to estimate 
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the IR luminosities of DOGs have primarily relied upon 
spec troscopic redshift s and the observed 24/xm flux den- 
sity. iDev et all (pOOSi ) use an empirical relation between 
the rest-frame 8/zm luminosity (computed from the ob- 
ser ved 24/xm flux density) and the IR luminosity, derived 
bv lCaputi et al.l ()2007^ . However, there is evidence from 
sources with F24^m > 0.25 mJy that methods based on 
only the 24/im flux densit y can overestimate the IR lumi- 
nosity by factors of 2-10 (jPapovich et al.ll2007D . Results 
from deep 70/im and 160/im imaging of a sub-sample 
of 24/im-bright DOGs are consistent with this, favoring 
hot-dust dominated spectral en ergy distribution (SED) 
templates like that of Mrk 231 (jTvler et al.ll2009[ ) which 
lead to estimates of the IR luminosity that are on the low 
end of the range in Lir/Ls conversion factors adopted in 
IDev et all mm . 

In this paper, we present 350/im and 1mm photom- 
etry of a sample of DOGs whose mid-IR spectral fea- 
tures (silicate absorption, power-law SEDs) suggest the 
presence of a strong active galactic nucleus (AGN). The 
primary goals of this study are to measure their IR lumi- 
nosities and constrain their dust properties, in particular 
the dust masses and temperatures. We also estimate stel- 
lar masses for the sources in the sample using published 
Hubble Space Telescope {HST) data and Spitzer lidraRed 
Array Camera (IRAC) catalogs fro m the Spitzer Dee p 
Wide-Field Survey (SDWFS; see lAshbv et al.l 120091) . 
Comparison of the stellar and dust masses potentially 
allows us to place constraints on the evolutionary status 
of these sources. 

In section [21 we present the details of the observations. 
Section [3] presents the DOG SEDs from 0.4/im to 1mm 
and IR luminosity measurements, constraints on the dust 
temperature, and dust and stellar mass estimates. In sec- 
tionlH we compare our results with similar studies of sub- 
millimeter galaxies (SMGs) and 5piteer-selected sources 
from the eXtragalactic First Look Survey (XFLS) and 
Spitzer Wide InfraRed Extragalactic (SWIRE) survey. 
We present our conclusions in section [5] 

Throughout this paper we assume a cosmology where 
Ho =70 km s-i Mpc~^, = 0.3, and Qx ^ 0.7. 

2. OBSERVATIONS 

2.1. Sample Selection 

IDev et al.l (poM ) identified 2603 DOGs in the 8.6 deg^ 
NDWFS Bootes field, selecting all 24/im sources satis- 
fying R - [24] > 14 (Vega mag) and ^24^1,1 > 0.3 mJy. 
We identified 12 DOGs with spectroscopic redshifts for 
follow-up 350/Ltm imaging (see Figure [T]) with the second- 
generation Submillimeter High Angular Resolution Cam- 
era (SHARC-II) at the Caltech Sub-millimeter Observa- 
tory (CSO). These targets were selected to have bright 
24/im flux densities (-F24^im ^ 1 rnJy) and a power-law 
dominated mid-IR SED (based on Spitzer/IKAC and 
24/im MIPS pho tometry; for details, see section 3.1.2 in 
IDev et all I2008D . Using the dee per IRAC observations 
from SDWFS (|Ashbv et al.l l2009). the fraction of DOGs 
qualifying as power-law sources ranges from « 10% at 
-^24^1,1 — 0.3 mJy to w 60% at F24^ni — 1 niJy. As 
shown in Figure [1] our sample spans a broad range in 
R - [24] color (« 14.5 - 17.5). Spitzer/IRS spectroscopic 
redshifts have been obtained for these sources based on 
the 9.7//m silicate absorption feature. Power-law con- 
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Fig. 1. — R — [24] color vs. 24/im magnitude for DOGs in 
the NDWFS Bootes field. Bottom and top abscissae show the 
24/im magnitude and flux density, respectively. Left and right 
ordinates show the color in magnitudes and the -F24(jm/i^o.7fim flux 
density ratio, respectively. Black dots and upward arrows show 
the full sample of DOGs, with and without an ij-band detection, 
respectively. The subsample studied in this paper is represented by 
red circles (open symbols show sources undetected in the _R-band 
data). Two sources observed by CARMA at 1mm are highlighted 
by a black square. This plot demonstrates that the sample studied 
in this paper probes the 24/im-bright DOGs over a wide range of 
R - [24] colors. 



tinua and silicate ab sorption are typical feattires of AGN 



dominated systems (iDonley et al.ll2007t Weedman et al.l 
l2006aHPonetta et al.ll2008t iBrand eraLll2008D . We note 
that such systems often exhibit intense star-formation 
concurrent with the g rowth of a super-massive central 
black hole (e.g.. IWang et al.ll2008f ). 

Details of our observations are presented in Table [TJ 
The effective integration time (column 8 of table [T]) 
represents the time necessary to reach the same noise 
level given a com pletely transparent atmosphere (see 
ICoppin et aT1l2008l for details). 

We observed two of the twelve DOGs at 1mm using 
the Combined Array for Research in ]VIillimeter-wave 
Astronomy (CARMA) interferometer to search for ther- 
mal emission from cold dust particles. These were pri- 
marily selected to have robust 350/im detections to en- 
hance the probability of detection at 1mm and, in the 
event of a non-detection, allow useful constraints to be 
placed on the dust properties. The two targets observed 
with CARMA are SST24 J142827.2+354127 (S2) and 
SST24 J143001. 9-^334538 (S3). Table[2]presents the date 
and integration time of the CARMA observations. 

2.2. SHARC-H SSOfim Imaging and Photometry 

The SHARC-II observations of the 12 target DOGs 
were carried out over the course of five separate observ- 
ing runs from 2005 January to 2007 May. Data were col- 
lected only when the atmospheric opacity was low and 
conditions were stable (r225GHz < 0.06). Pointing, focus 
checks, and calibration were performed every hour, using 
ULIRG Arp 220 as a calibrator (Faso^m = 10.2 ± 1.0 Jy). 
Other secondary calibrators (CIT6, CRL2688, 3C345) 
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TABLE 1 
SHARC-II 350/iM Observations 



Source Name 


ID 


RA 
(J2000) 


DEC 
(J2000) 


z 


UT Year-Month 


(hr) 


*cff 

(min) 


(mJy) 


SST24 J142648. 9+332927 


SI 


14:26:48.970 


+33:29:27.56 


2.00^= 


2006-Apr 


1.1 


4.0 


20 


SST24 J142827.2+354127 


S2 


14:28:27.190 


+35:41:27.71 


1.293=* 


2005-Apr /2006-Apr/2007-May 


3.3 


11.9 


10 


SST24 J143001. 9+334538 


S3 


14:30:01.910 


+33:45:38.54 


2.46° 


2005- Apr/2006- Apr/2006-May 


3.6 


11.4 


8 


SST24 J143025. 7+342957 


S4 


14:30:25.764 


+34:29:57.29 


2.545* 


2006-Apr 


1.0 


2.0 


26 


SST24 J143135. 2+325456 


S5 


14:31:35.309 


+32:54:56.84 


1.48== 


2007-May 


1.2 


2.0 


30 


SST24 J143325. 8+333736 


S6 


14:33:25.884 


+33:37:36.90 


1.90== 


2006-Apr 


0.3 


1.2 


28 


SST24 J143411. 0+331733 


S7 


14:34:10.980 


+33:17:32.70 


2.656S 


2005-Jan 


3.8 


12.7 


8 


SST24 J143447.7+330230 


S8 


14:34:47.762 


+33:02:30.46 


1.78== 


2007-May 


1.8 


10.5 


10 


SST24 J143508. 4+334739 


S9 


14:35:08.518 


+33:47:39.44 


2.10'= 


2006-Apr 


1.0 


0.5 


35 


SST24 J143539. 3+334159 


SIO 


14:35:39.364 


+33:41:59.13 


2.62<= 


2005- Apr/2006- Apr 


2.5 


3.1 


15 


SST24 J143545. 1+342831 


Sll 


14:35:45.137 


+34:28:31.42 


2.50== 


2006-Apr 


0.8 


3.2 


22 


SST24 J143644.2+350627 


S12 


14:36:44.269 


+35:06:27.12 


1.95== 


2006-Apr/2006-May 


2.1 


6.7 


10 



Actual on-sourcc integration tim 
Effective integration time for a t 
' Redshift from S'pitzer/IRS (Higdo 
* Redshift from Keck DEIMOS [H] 
' Redshift from Spit^er/IRS 

Redshift froi 
' Redshift fro 



et al. 



ent atmosphe: 
in prep). 



k et al.irWTBI ). 

k DEIMOS (Dey et al., in prep.). 
:k LRIS jl3ev et am^OOSI) . 



TABLE 2 
CARMA iMM Observations 



ID 


UT Year-Month 


tint (hr)=' 


S2 


2008-April/May 


10.3 


S3 


2008-April/May 


7.5 



On-sourcc integration time. 



were used occasionally to verify the Arp 220 calibrations. 
A non-connecting Lissajous pattern was used to modu- 
late the telescope pointing with amplitudes of 15"-20" 
and periods of 10-20 s. The observations made use of the 
CSO Dish Surface Optimization System to optimize the 
dish surface accuracy and beam efficiency (jLeong et al.l 
[2001 . 

Data were reduced using the Comprehensive Reduc- 
tion Utility for SHARC-II (CRUSH) software package 
with the 'deep' option to optimize the signal-to- noise ra- 
tio (S /N) for faint (< 100 mJy) point sources (jKovacsl 
I2006f ). The output map has a pixel scale of 1'.'62 pixel" 
and is smoothed with a 9" gaussian beam, resulting in 
an effective image FWHM of 12'.'4. 

A 20" diameter aperture was used for photometry to 
compute the instrumental flux density of each source. 
The sky level and photometric uncertainty were com- 
puted by measuring the mean and RMS in «10 off-source 
20" diameter apertures. The same procedure was ap- 
plied to the calibration images, and a scaling factor was 
derived that converts the instrumental flux density to a 
physical flux density (using this method, no subsequent 
aperture correction is required as long as both the science 
and calibration targets are unresolved and measured in 
the same aperture). 

The aperture photometry is consistent with peak 
flux density measurements in all but one source, 
SST24 J142648.9+332927. This source has a radial pro- 
file that is significantly more extended than the point 
spread function of the final map, which results in the 



peak flux underestimating the aperture flux density mea- 
surement. The extended structure in the image is more 
likely to be noise than signal, so in this case we report 
the peak flux measurement, which is formally a non- 
detection. 

Flux boosting of low S/N sources can be an impor- 
tant effect in wide-field surveys where source positions 
are not known a priori (e.g., Coppin et al.' 2005|). How- 
ever, because we know our source positions at the < 1" 
level (from MIPS and IRAC centroids), flux boosting is 
not a significant effect, and so we do not apply any such 
corrections to o ur measuremen t s. Ou r appro ach follows 
that a dopted bv lLaurent et al.l ()2006D and K ovacs et all 
(|2006D in their SSO^im follow -up imaging of SMGs. 

2.3. CARMA 1mm Imaging and Photometry 

The CARMA observations were obtained between 2008 
April 7 and May 1 in the C-array configuration (beam- 
size is « 2 X 1 sq. arcsec). A total of 7.5 hours and 10.3 
hours of integration time in good 1mm weather condi- 
tions were spent on sources SST24 J143001. 9+334538 
(S3) and SST24 J142827.2+354127 (S2), respectively. 
These sources were selected primarily because of their 
robust (S/N> 4.5) detections at 350/im. In addition, 
source S3 is detected wit h Spitzer/MlPS at TOfim and 
160^m (jTvler et al.ll2009[ ). w hile S2 is the sub ject of a 
detailed spectroscopic study (|Desai "eFaIl[2006h . 

System temperatures were in the range 250-400 K. A 
correlator configuration was used with three adjacent 
15x31 MHz bands centered on 220 GHz, the frequency 
at which the CARMA 1mm receivers are most sensitive. 
The quasar J1310+323 (chosen for its spatial proximity) 
was observed every 15 minutes for amplitude and phase 
calibration. Quasars 3C 273 and MWC 349 were used 
for pointing, pass-band calibration, and flux calibration. 

Data were reduced using the Multichannel Image Re- 
construction, Im age Analysis, an d Display (MIRIAD) 
software package (jSault et al.lll995D . Visual inspection of 
visibilities as a function of baseline length allowed us to 
identify and flag spurious data. A cleaned map was gen- 
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erated for each track of integration time (ranging from 
1 to 5 hours in length) and these tracks were coadded 
together using a weighted mean to obtain a final im- 
age of these sources. No detections were found in either 
case. Both sources are unresolved in the IRAC images, 
and NICMOS/F160W imaging of S3 indicates a size of 
< 0'.'5, implying that it is very unlikely any emission was 
resolved out by the interferometer. The imstat routine 
from MIRIAD was used to determine the noise level in 
the co-added images where we expected to see the source. 
TableOshows the photometry from 24/im to 1mm, where 
available. Non-detections are given as 3-tT upper limits. 

2.4. Optical, near-IR, mid-IR, and far-IR Photometry 

The optical and near-IR photometry used in this pa- 
per to estimate stellar masses are based on high spatial- 
resolution HST data (WFPC2 /F606W, ACS/F814^y , 
and NIC2/F160W) published in lBussmann et all (|2009D . 
The HST data allow the separation of an unresolved nu- 
clear component (flux on scales < 1 kpc) from a more 
spatially extended component. Because the AGN con- 
tribution in the rest-frame UV to optical is uncertain, 
the photometry of the extended component is used here 
(measured with 2" diameter apertures) to ensure that our 
stellar mass estimates are not biased by the presence of 
an obscured AGN (for de tails on how the extend ed com- 
ponent is computed, see iBussmann et al.l 120091 ) . Addi- 
tionally, 4" diameter aperture photometry in the optical 
{Bw, R, and /) from the NDWFS is shown in the SEDs 
of the objects in this sample (details on h ow the pho- 
tome try is computed may be found in Bus smann et al.l 
[20091) . 

The mid-IR photometry used in this paper are from 
the publicly available Data Release 1.1 (DRl.l) cata- 
logs from the SDWFS IRAC coverage of the Bootes 
field (Ashby et al. 2009). The SDWFS catalogs incor- 
porate the earlier IRAC Shallow Survey of the Bootes 
field undertaken by th e IRAC guaranteed tim e observa- 
tion (GTO) programs (jEisenhardt et al.ll2004D . We iden- 
tified IRAC counterparts of the DOGs in this paper from 
the SDWFS catalogs using a 3" search radius centered 
on the 24/im position (the MIPS 24/im 1-a positional 
uncertainty is 1'.'2). All of the DOGs in this paper have 
IRAC counterparts, detected at > 4tT in all four IRAC 
channels. We use the 4" (rather than the 6") diameter 
aperture photometry from the DRl.l SDWFS catalog 
to reduce contamination from nearby sources. We note 
that aperture corrections derived from isolated, bright 
stars have been applied to the SDWFS catalogs. 

Finally, 24, 70, and 160/im data over 8.61 deg^ of the 
Bootes field are available from GTO programs. The data 
were reduced by the MIPS GTO team and reach la rms 
depths of 51 /tJy, 5 mJy, and 18 mJy at 24, 70, and 
160/fm, respectively. Details of the GTO surveys, such 
as mapping strategy, data reduction, and source cata- 
logs, will be discussed elsewhere. In addition, several of 
the DOGs in this paper were targeted for deeper MIPS 
photometric observations by Spitzer General Observer 
program 20 303 (P.I. E . LeFlo c'h), and the results are 
report ed in iTvler et all (|2009D . We use the iTvler etld] 
(|2009f ) measurements where they are available. 

3. RESULTS 



In this section, we present SEDs from 0.4/tm to 1mm 
for each source in our sample and compare with local 
starburst (M82) and ULIRG (Arp 220 and Mrk 231) tem- 
plates. Our approach is to artificially redshift the local 
galaxy templates and normalize them to match the DOG 
photometry at observed-frame 24/tm. This allows a sim- 
ple, qualitative comparison of DOGs and galaxies with 
properties ranging from warm dust, star-formation dom- 
inated (M82), to cool dust, star-formation dominated 
(Arp 220), to warm dust, AGN-dominated (Mrk 231). 
We will use the SED that provides the best fit over the 
sampled wavelength range (observed optical through sub- 
mni) to estimate the IR luminosities of the DOGs. 

Later in this section, we use our limits at 1mm from 
CARMA to place constraints on the dust temperatures 
and limits on the dust masses of DOGs. Finally, we use 
HST and Spitzer/IRAC data to estimate stellar masses 
of DOGs. 

3.1. Qualitative SED Comparison 

Figures [2] and [3] show the SEDs of the sources, divided 
respectively into those with and those without detections 
at 350/tm. Note that the rest-frame UV photometry 
for source S7 is contaminated by emission from nearby 
sources and will be treated in more detail in a future 
paper (Prescott et al., in prep.). 

Overplotted in each panel are M82 CSil va et al.lll998l )'^. 
Mrk 231 (Chary 2008, p rivate communication), and 
Arp 220 (|Rieke et al.ll2008D templates, placed at the ap- 
propriate redshift and scaled to match the flux density 
observed in the MIPS 24/tm band. The scaling factors 
derived for the three templates range over 200-900, and 2- 
10, 70-700, respectively (the deep silicate absorption fea- 
ture in Arp 220 and the strong PAH emission feature of 
M82 make the scaling factors closer to each other than a 
simple estimate based on the ratio of the IR luminosities 
would imply). For the Arp 220 and Mrk 231 templates, 
we have interpolated the spectrum in the UV to match 
Galaxy Evolution Explorer photometry (in the case of 
Arp 220) and International Ultraviolet Explorer data 
as well as HST Faint Object Spectrograph data (in the 
case o f Mrk 231: lHutchings fc NefjllgSTl : [Gallagher et all 
[200l . 

These templates were chosen because they sample a 
range of dust temperatures and AGN/star burst contri- 
butions. M82 is one of the closest (dL — 3.86 Mpc) 
galaxies undergoing a starburst, as it was triggered by 
a recent interaction with M81. Although it is less lumi- 
nous than DOGs (Lir « 6 x 10^° Lg; ^ iSanders et all 
I2003f ). its nucleus is dominated by a war m dust compo- 
nent (Tdust = 48 K; IHughes et a"l]|1994D . Arp 220 is a 
nearby (rff, = 77.3 M pc) ULIRG (Lir « 1.6 x lO^^ L©; 
ISanders et al.l 120031) dominate d by cold dust (Tdust = 
35 K: IRigopoulou et all Il996f) . Mrk 231 is another 
nearby (dr, = 175.1 Mpc) ULIRG (Lir « 3.2 x lO^^ L©; 
Sanders et al. 2003), but has a warm dust (Tdust = 51 K; 
I Yang fc Phillips ,2007.) SED dominated by an obscured 
AGN. 

Qualitatively, the Mrk 231 template provides a much 
better fit than the Arp 220 template to the 350/tm pho- 

We use a slightly updated SED obt ained from 
|http://adlibitum. oat. ts. astro. it/silva/grasil/modlib/m odlib.html| 
° -Lib, is the luminosity integrated over 8-1000/xm 
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TABLE 3 
Photometry'' 





R - [24] 


-^24^111 




-F'leo/jm 




-F'l.2mm 




Source Name 


(Vega mag) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


SI 


>16.1 


2.33±0.07 






< 66 






S2 


15.7 


10.55±0.13 




<45': 


74±13 


<1.5 




S3 


>17.4 


3.84±0.06 


9.3±2.3<i 


65±lld 


41±13 


<1.8 


0.42 ± 0.04 


S4 


15.1 


2.47±0.05 






<81 




0.23 ±0.03 


S5 


14.5 


1.51±0.05 






<100 




0.54 ±0.12 


S6 


15.4 


1.87±0.06 






<137 




0.20 ±0.03 


S7 


12.4" 


0.86±0.05 


<25'' 


<90'' 


37±13 






S8 


>16.8 


1.71±0.04 






45±12 
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Photometry from lTvler et al.l Il2009f) . 

^ i^-band photometry includes diffuse emission from Ly-a nebula jPev et al.ll200'5l) . 
f Photometry from lDev et al.l 1 12005ft . 
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Fig. 2.— SEDs of 5 DOGs detected by SHARC-II at 350Atm. Dotted, dashed, and dot-dashed lines show the Mrk 231, Arp 220, and M82 
template SEDs, respectively, placed at the appropriate redshift and scaled to match the observed 24/im flux density. Red horizontal lines 
show the 5-0- sensitivity limits (ignoring confusion) from the planned wide-field Herschel surveys at 250, 350, and 500/im, and SCUBA-2 
surveys at 850/im. The cool dust BED of Arp 220 significantly overpredicts the SbOfira fiux density in all cases. The warm dust SED of 
M82 provides a better fit in the far-lR, but Mrk 231 provides the best fit in both the far-IR and the optical. 
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Fig. 3. — Same as Figure [2] except showing SEDs of 7 DOGs not detected by SHARC-II at 350/im. The limits at 350/im are all 
inconsistent with the Arp 220 SED. The 250/im channel of Herschel should be very efficient for detecting power-law DOGs in wide-field 
surveys, assuming that Mrk 231 is an appropriate representation of the far-lR SED. 
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tometry in every case. M82 fits the 24/xm and 350/L(m 
photometry reasonably well (although not as well as 
Mrk 231), but it fares poorly in the mid-IR and op- 
tical, where a strong stellar component in M82 is not 
seen in the DOGs in this sample (which are dominated 
by a power-law in the mid-IR) . Additionally, M82 shows 
strong PAH emission which is not seen in the power-law 
DOGs. 

The red horizontal bars in Figures [2] and [3] show 5- 
a limits (ignoring confusion) from planned wide-field 
(> 8 deg2) surveys with the Herschel Space Observa- 
tory (shown for the channels at 250, 350, and 500/im) 
and with the Sub-mm Common-User Bolometer Array-2 
(SCUBA-2) instrument at 850 /xm. Most of the power- 
law DOGs studied in this paper have SEDs that peak 
around observed-frame 250/im, which is where the Her- 
schel wide-field maps will be the deepest (5-(t limit of 
14 mJy). If all of the 24/xm-bright DOGs are detected 
at 250/um in the two wide-field surveys that are planned 
to reach the depths assumed here (Lockman Hole east, 
11 deg^; Extended Chandra Deep Field South, 8 deg^), 
then a total of « 600 power-law DOGs should be de- 
tected in the 250/im Herschel catalogs of these two fields. 
The SCUBA-2 surveys of these fields should be deep 
enough (5-tT limit of 3.5 mJy at 850/im) to detect many 
of these sources, allowing dust temperature constraints 
to be placed on a statistically significant sample of these 
rare, important objects. 

Figure m shows all of the DOG SEDs on the same plot, 
normalized by the rest-frame 8/im flux density (which 
is estimated from the observed-frame 24/im flux den- 
sity by assuming a power law of the form F^, oc i>°', 
where a = -2). The SEDs M82, Arp 220, Mrk 231, 
and a composite SMG template SED spanning mid-IR 
to sub-mm wavelengths are also shown. The compos- 
ite SMG template is derived from bright (i^ssoAmi > 
5mJy) SMGs from the Great Observatories Origins Deep 
Survey North (GO ODS-N) field with mid-IR spectra 
(|Pope et a l. 2008bl. 

One striking feature of this plot is the steep slope 
shown by DOGs in the rest-frame 1-4/im. Whereas 
Mrk 231, M82, and Arp 220 all exhibit a bump in the 
l-2fim regime, no such feature is apparent in the DOG 
SEDs. This could be due to the presence of an obscured 
AGN outshining the stellar light, in the rest-frame near- 
IR. This is in constrast to rest-frame UV and optical 
wavelengths, where HST ima ging has revealed that stel- 
lar light appears to dominate (jBussmann 
noted previously, DOGs have far-IR to mid-IR flux den- 
sity ratios more similar to Mrk 231 than Arp 220. The 
composite SMG template overpredicts the far-IR flux for 
a given mid-IR flux in all cases where we have 350/im 
detections. We note that adding an additional warm 
dust component (Tdust = 350 K; possibly powered by 
an AGN) to the c omposite SMG SED (as was done in 
iPope et all l2008a[ ) improves the quality of the fit over 
the rest-frame 8-100/^m. However, this composite SMG 
-I- AGN template retains a strong cool dust (Tdust ~ 30K) 
component that over-predicts the amount of emission at 
1mm. If this type of SED was appropriate for the power- 
law DOGs investigated in this paper, they would have 
been easily detected by CARMA. 

An alternative way of displaying this information is 



shown in Figure [51 In each panel, the flux density ratio 
far-IR:mid-IR is plotted as a function of the flux density 
ratio mid-IR:optical (F24^im/-Fo.7^m)- The top two panels 
show -F'350/^m/-F24;^m On the y-axis, while the bottom two 
panels show Fi200fj.m/ F24:tim on the y-axis. SMGs and 
various Spitzer-selected sources are shown in the plots, 
divided into those that are detected in the (sub-)mm on 
the left and those that are not detected on the right. 
The SMG XFLS , and SWIRE R-h and data come from 
'Dve et_alj (|2008| ) . lYan et all ()2007f ). and lLonsdale et al.l 
(2009), respectively. For SMGs without detections at 
1200/im, Fi2oo;im is estimated usin g the 850/im flux den - 
sity and the dust temperature from lCoppin et all (|2008f) . 
and is represented by a red cross symbol. Dotted, dashed, 
and dot-dashed lines indicate the evolution of Mrk 231, 
Arp 220, and M82, respectively, on this diagram over red- 
shifts of 1 — 3. Compared to SMGs, DOGs in this sample 
have redder flux density ratios in the mid-IR:optical but 
bluer far-IR:mid-IR ratios. This cannot be explained by 
an enhancement of the 24/im flux due to PAH emission, 
since mid-IR spectra of these DOGs show power-law con- 
tinua with silicate absorptio n and weak or absent PAH 
emission features (iHouck et a l. 2005). Instead, the most 
likely explanation is that ob scured AGN emission bo osts 
the mid-IR continuum (e.g.. lRieke fc Lebofskvl[T98lf) rel- 
ative to both the optical and far-IR. 

3.2. IR Luminosities 

In this section, we provide the best available estimates 
of the total IR luminosities (Xir; 8-1000/im rest-frame) 
of the sample in this paper based on the 350/im imaging. 
We then compare these with estimates based solely on 
the 24/im flux density, and also with estimates of the far- 
IR luminosity (Tfir; 40-500/im, rest-frame) based on a 
modified black-body which has been scaled to match the 
sub-mm photometry. 

The qualitative SED comparison from section [33] sug- 
gests that Mrk 231 provides a reasonable flt to the far- 
IR photometry. In addition, analysis of 70/xm and 160/im 
photometry of a sample of these types of AGN-dominated 
DOGs has suggested that Mrk 231 provides a reason- 
able approximation of the full SED (j Tyler et al.l I2009I . 
see table [3] for overlap between that study and this 
one). Therefore, as the best measure of Lir, we inte- 
grate (over 8-1000/xm rest-frame) a redshifted Mrk 231 
template which has been scaled to match the observed 
350/im flux density (or 3-a limit, in the case of a non- 
detection). These values are tabulated in the first column 
of Table [H For sources with detections at 350/im, Lir is 
in the range of (2.0 — 2.6) x 10"'^^ Lq, with a median value 
of 2.2 X 10^3 Lq. 

The second column of Table [4| shows the rest-frame 
8/fm luminosity, j^Li,(8/im), while the third column 
shows the total IR luminosity (8-1000/fm, rest-frame) 
based on the spectroscopic redshift and an empirically 
determined relationship between i/L ^{8^m) and Lir: 
Lir = 1.91L^ °^ ("Caputi et all l2007^. This approach was 
used bv lDev e~ al. ( 2008) in determining the contribution 
of DOGs to the total IR luminosity density of all z ~ 2 
galaxies. Lir values range over (0.5-8.1) x 10^^ Lq, with 
a median value of 2.3 x 10^^ Lq. The Caputi et al. de- 
rived value for Lir is consistent with the 350/im-based 
estimate (or 3-a limit, in the case of non-detections) in 
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Fig. 4. — Optical through sub-mm SEDs of DOGs in the SHARC-II sample. Flux densities have been normalized by the rest-frame 8/im 
flux density, computed from the observed 24^m flux density. Of the three local galaxy templates shown, Mrk 231 provides the best fit over 
the rest-frame UV through sub-mm range because it has a warm dust SED (unlike Arp 220) and because it lacks a strong stellar component 
(unlike both Arp 220 and M82). None of the template SEDs match the steepness of the rest-frame near-IR photometry of the power-law 
DOGs. This may indicate that obscured AGN dominate stellar emission to a greater extent in power-law DOGs than in Mrk 231. 



6/12 of the power-law DOGs studied here. In the re- 
maining half of the sample, the Caputi et al. relation 
overestimates Lir in 5/6 targets. In only one DOG (S8) 
is the 350^m emission brighter than would be expected 
based on the 24^m flux density, redshift, and the Ca- 
puti et al. relation. This implies that measurements of 
the IR luminosity density of DOGs relying solely on the 
24/im flux density will tend to overestimate their true 
contribution, consistent with what has been found in a 
recent stud y of faint (J24^/,m ^ 100 — 500/iJy) DOGs in 
GOODS-N (|Pope et al.ir2008aD . Quantifying the extent 
of this effect will require much larger samples of DOGs 
with sub-mm measurements, the kind that will result 
from wide field surveys with Herschel and SCUBA-2. 

Finally, the last column of Table |4] shows FIR lumi- 
nosities computed from the integral over 40-500/im (rest- 
frame) of the best-fit modified black-body (described in 
more detail in section [33)1 . These values are tabulated 
only for those sources with CARMA 1mm imaging. We 
find LpiR values of « 10^"^ Lq, implying LiR/ipiR ~ 3. 
In contrast, Mrk 231 has Xir/Lfir, ~ 2, underscoring the 
fact that the IR luminosity of these DOGs is dominated 
by mid-IR emission rather than FIR emission. 

3.3. Constraints on Dust Properties 

Additional constraints can be placed on the nature of 
the cold dust emission from the two sources (S2 and 
S3) with CARMA. The 1mm non-detections imply warm 
dust temperatures. If we compute the predicted flux den- 
sity at 1mm based on the observed 24/Ltm flux density and 



assuming the three local galaxy SED templates (M82, 
Arp 220, and Mrk 231) used in the previous sections, we 
find values of 4.4, 740, 1.8 mJy, and 5.8, 120, 2.7 mJy 
for the two sources respectively. For 82, this implies that 
our 3-tT limit (< 1.5 mJy) is close to the level we would 
expect if Mrk 231 is an appropriate SED, while the other 
two SEDs are clearly inconsistent with the data. For S3, 
the limit (< 1.8 mJy) is inconsistent with each of the 
SEDs, implying a warmer dust temperature than even 
Mrk 231. 

A more quantitative approach is to minimize the resid- 
uals between the 160/xm, 350/.tm, and 1mm data and that 
expected from a modified black body. Doing this, we can 
constrain the dust temperature, Tdust, and the dust emis- 
sivity index, f3, of each target. The best-fit quantities and 
their uncertainties are estimated using a bo otstrap tech- 
nique that mimics the procedure used by iDunne et al.l 
(2000;). Briefly, for each flux density measurement, a set 
of 100 artificial flux densities are generated using a Gaus- 
sian random number generator. The mean value and dis- 
persion of the distribution of artificial flux densities are 
set by the measurement value and its la uncertainty, re- 
spectively (for non-detections, we assume a mean value 
of and force artifical flux densities to be positive). Each 
artificial SED is used to construct a distribution of best- 
fit scaling factors and associated xt values for a modified 
black-body with a given combination of f3 and Tdust • 

Figure [6] shows the median xt contours for the grid 
of Tdust and (3 values we have sampled in the model 
fitting for each source observed by CARMA. The con- 
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Fig. 5. — Top and bottom panels show 350fim/24fim and 1200/im/24/im flux density ratios, respectively, as a function of R-[24] (Vega 
mag) color. For clarity, sources are separated into those detected at 350/im or 1200^m on the left and those that are not detected on 
the right. Objects that qualify as a DOG (R - [24] > 14) are shown with a filled symbol. Black circles indicate DOGs in Bootes. Red 
squares show measurements of SMGs l lCoppin et al.l [2008!). while red crosses show predicted valu es based on 850/jm photometry (see text 
for details). Green stars show Spitzer-selected bump sources from SWIRE (jLonsdale et al.?'200(3). Grange triangles (teal triangles) show 
similarly identified sources from the XFLS dom inated by PAH emission (silicate absorption) features (Sajina et al. 20Q3|)- Magenta inverted 
triangles show XFLS sources from iSaiina et al.l | |200 8). Finally, dotted, dashed, and dot-dashod lines show the evolution of Mrk 231, Arp 220, 
and M82 in this parameter space over redshift 1 to 3. The DOGs studied in this paper have some of the reddest R — [24] colors and lowest 
far-IR/mid-IR flux density ratios of other z ~ 2 ULIRGs like SMGs or other Spifaer-selected sources. 



TABLE 4 
Luminosities 



Source Name 


(1012 Lq) 


(1012 Lq) 


(1012 Lq) 


(1012 Lq) 


SI 


< 32<= 


2.2 


23 




S2 


26 ±5^1 


2.4 


25 


9.1 


S3 


22 ±7^1 


7.2 


81 


10 


S4 


< 44<= 


5.2 


57 




S5 


< 40<= 


53 


5.1 




S6 


< 64<= 


1.5 


15 




S7 


21 ± 


2.1 


22 




SB 


20 its'! 


1.1 


11 




S9 


< 75° 


2.9 


31 




SIO 


< 28'= 


6.2 


69 




Sll 


< 33<= 


3.9 


42 




S12 


< 16<= 


2.0 


21 





^ Integral over 8-1000/xm of rcdshifted Mrk 231 template normalized at 

350/^m. 

^ Estimated from lyL^ {Sfim) - Li^ relation from fUaputi et al.| j2007D . 

Integral over 40-500/im of best-fit modified black-body (only sources 
with CARMA 1mm data). 

Uncertainties shown reflect 350/xm photometric uncertainties. Ad- 
dtional systematic uncertainties associated with the adoption of a 
Mrk 231 template are not included. 

3(7 upper limits. 
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Fig. 6. — x2 contours based on modified black-body fits to sources 
with CARMA data. Lines indicate 80%, 95%, and 99.9% confi- 
dence interval s. Also shown are Tjust distributions for SMGs with 
350Atm data I IKovacs et al.|[2003 : [Coppin et all 12001) . The 95% 
confidence levels for S2 and S3 suggest dust temperature limits 
that would place them in the warmest 50% and 15%, respectively, 
of SMGs. 

tours show the degeneracy between (3 and Tdust, and il- 
lustrate why a perfect fit is not possible in spite of the 
fact that a model with three parameters is being fit to 
three data points (i.e., the input parameters are not fully 
independent). The distributio n of Tdust foun d for S MGs 
based on SSO /xm imaging from lCoppin et all (j2008f ) and 
IKovacs et al.l ([2006) is displayed in the lower panel. Two- 
sided Kolmogorov-Smirnov tests suggest that the distri- 
bution of Tdust (assuming /3 = 1.5) for each CARMA 
target is highly unlikely to be drawn from the same par- 
ent distribution as the combined sample of SMGs (4% 
and 0.008%, for the two CARMA targets respectively). 
In both of the CARMA targets, warmer dust temper- 
atures are required due to the non-detection at 1mm. 
The data cannot rule out models with higher values of 
/3, but models of dust grains as well as Galactic and 
extra-galactic observations consistently suggest (3 = 1 — 2 
(jHildebrand 1983; Dunne & Ealcs 2001). For (i = 1.5, 
we reject Tdust < 33K and 45K at the 95% confidence 
level, in S2 and S3 respectively. For /? = 2, we can reject 
Tdust < 25 K and 37 K. 

3.4. Dust Masses 

Assuming optically thin sub-mm emission, cold dust 
masses can be estim ated from the 350/xm photometry 
(jHughes et al.lll997f) : 



Md„ 



SohsD^ 



1 + z Kf''*^B{v"'''\ 



Tdust) 



(1) 



where S'obs is the observed 350^m flux density, and n^^^^ 
and 5(1^''°'^*, Tdust) are, respectively, the values of the 
mass absorption coefficient and black-body function at 



the rest frequency 1/^°^^ and dust temperature Tdust- The 
approp riate value for Kd is uncertain to at least a factor 
of two (jDunne et al.ll200j); we use a k'^J^^ value interpo- 



lated from Drainc (2003) {<k 



20 cm^ g- 



The results from section 13.31 suggest that Tdust > 35 — 
60 K for two of the sources. Adopting the average of 
these limits (Tdust — 45 K) for all of the sources in this 
sample, then dust mass limits are in the range (4.1 — 
7.9) X 10® Mq (median value of 5. 1x10® M0) for the five 
objects with detections at 350/im. The 3(t upper limits 
on the dust masses of the remaining sample range from 
(3 — 15) X 10® Mq. Warmer values of Tdust would lead 
to smaller inferred dust masses (e.g., increasing the dust 
temperature by 10 K implies «50% lower dust masses). 

The dust masses are presented in TableJS] 

The se values agree with those of iBussmann et al.l 
(|2009[ ). where dust masses of a sample of 31 power- 
law dominated DOGs were estimated using predicted 
850^m flux densities based on Mrk 231 templates and 
the measured 24/xm flux density. In previous work, 
we assumed a Tdust of 75 K and found a median dust 
mass of 1.6x10® Mq. This is consistent with the notion 
that Mrk 231 accurately characterizes the far-IR SED 
of power-law DOG s, as described in section 13.11 and in 
iTvler et al.l (p009h . 

Finally, assuming a gas mass to dust mass ratio of 120 
(as was fou nd in a study of the nuclear regions of nearby 
LIRGs; see lWilson et al.. 2008). then the gas masses can 
be estimated as well. Using the assumed ratio, we find 
gas masses of (5 — 10) x 10^° Mq (median value of 
6x10"^° Mq) for the detected objects and gas mass 3a 
limits of (4—18) x 10^^ Mq in the remaini ng samp l e. We 
caution that this is very uncertain; Kovacs et al.l ()2006f ) 
report a gas mass to dust mass ratio of ~ 60 for SMGs, 
assuming k™''' = 15 cm^ g"-*^. If this gas to dust mass 
ratio is appropriate for our sample, then the implied gas 
masses will be a factor of two lower (2 — 5 x 10^" Mq) 

3.5. Stellar Masses 

In this section, we describe the methodology and 
present estimates for the stellar masses of the DOGs in 
this sample. 

3.5.1. Methodology 

To estimate stellar masses, we rely on Simple 
Stellar Population (SSP ) template SEDs from the 
iBruzual fc Charlotl ()2003f ) population synthesis library. 
All models used here have ages spaced logarithmically 
from 10 Myr up to 1 Gyr, solar metallicity, a Chabrier 
initial mass function (IMF) over the mass range 0.1 — 
100 Mq (Chabrier 200 3), and use th e Padova 1994 evo- 
lutio nary tracks (iGirardi et al.lll99^ . The reddening law 
from lCalzetti et al.| ( 2OO0 ) is used between 0.12 — 2.2 ^m 
and that of lOraind ( 200S ) for l onger wavelength s . Thi s 
method is similar to that used in lBussmann et aTl (120091 ). 

For sources at z ^ 2 whose mid-IR luminosity is dom- 
inated by stellar light, IRAC photometry samples the 
SED over the wavelength range where emission from 
asymptotic and red giant branch stars as well as low-mass 
main-sequence stars produces an emission peak at rest- 
frame 1.6/zm. In such cases, for given assumptions re- 
garding the star-formation history, metallicity, and IMF, 
stellar mass estimates can be obtained via stellar pop- 
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ulation synthesis modeling. One goal of this work is to 
estimate stellar masses using self-consistent modeling of 
photometry measured at similar wavelengths for a variety 
of z ^ 2 dusty galaxies. Therefore, we apply this method 
to determine stellar masses in SMGs as well as XFLS and 
SWIRE sources. The IRAC data for e ach of these galax y 
population s come s fro m, respectively, Dve et all (|2008D . 
iLacv et al.l (|2005D . and lLonsdale et al.l 1200^ . 

The DOGs studied in this paper have mid-IR SEDs 
that are dominated by a power-law component, suggest- 
ing that obscured AGN emission is overwhelming the 
stellar flux at these wavelengths. The shape of the mid- 
IR SED therefore provides limited constraints on the stel- 
lar population and additional information is needed to es- 
timate the stellar mass of these sources. To overcome this 
challenge, SSP models were fit to high spatial-resolution 
7757 photometry in the rest-frame UV (WFPC2/F606W 
or AC S/F814W) and r es t-fram e optical (NIC2/F160W) 
from iBussmann et al.l (|2009D . Two sources cur- 
rently lack HST data (SST24 J142827.2+354127 and 
SST24 J143411.0-f 331733^) and so are excluded from 
this analysis. 

In principle, rest-frame near-IR data offer a means to 
estimate the SSP age and Ay independently, since rest- 
frame optical and near-IR photometry sample the SED 
above the 4000 A break while the rest-frame UV photom- 
etry samples galaxy light below the 4000 A break. How- 
ever, these data come from the IRAC 3.6/im images of the 
Bootes field {Spitzer Deep Wide-Field Survey; Ashby et 
al., in prep.), where the spatial resolution is insufficient to 
resolve the nuclear source from the extended galaxy com- 
ponent. In this case, there are only limited constraints on 
the amount of non-stellar (i.e., obscured AGN) emission 
at 3.6/im. We have explored the effect of this uncertainty 
on the fitting process by artificially reducing the 3.6/im 
flux by 50% (corresponding to the situation where the 
3.6/im emission is equal parts starlight and AGN) and 
re-analyzing the data. Comparing results, we find that 
higher AGN fractions imply younger inferred SSP ages 
and higher Ay values. Because the AGN fraction in these 
sources is currently unknown at 3.6/im, we are unable to 
constrain both the age and Ay independently. 

Although the usage of photometry at different wave- 
lengths is not ideal for the purposes of comparing stellar 
masses between galaxy populations, this method remains 
valuable because the models being fit to the data are the 
same for each galaxy population. Indeed, recent work has 
suggested that the dominant source of systematic uncer- 
tainty in stellar mass estimates of K-selected galaxies at 
z ~ 2.3 is the use of different stellar population synthe- 
sis codes (jMuzzin et al.l[2009t ). and that these systemat- 
ics often dominate the formal random uncertainty. As 
long as the parameters of the model used here (such as 
the IMF, star-formation history, metallicity, etc.) do not 
vary from population to population, then the comparison 
presented here should be valid in a global sense. 

3.5.2. Stellar Mass Estimates 

Figure [7] shows xl contours for a grid of SSP ages and 
Ay values. The contours trace fines of 80%, 95%, and 
99.9% confidence intervals allowed by the photometric 

^ HST data exist for this source but will be presented in a 
separate paper (Prescott et al., in prep.) 



uncertainties (estimated using a bootstrap method sim- 
ilar to that outlined in section [373)1 . The solid grey lines 
trace iso-mass contours and show the range in stellar 
masses allowed by the photometric uncertainties. The 
best-fit SSP model parameters (Mgtar and x1) are printed 
in each panel and shown in Table El 

The stellar masses in the sample range from (1 — 20) x 
10^° Mq, with a median value of 2 x 10^° M©. The xl 
values range from 0.37 — 1.16, with a median value of 
0.69. In one case (SST24 J143001. 9+334538, or S3), the 
photometric uncertainty is so large that < 1 over the 
full range of Ay and SSP age that we have sampled and 
so the range of acceptable fits is very large. For this 
source, we quote the 3a upper limit on the stellar mass 
based on the photometric uncertainty. 

The ratio of the stellar to gas mass, C = Mstar/-^^gas, 
is a measure of the evolutionary state of the galaxy, with 
larger values indicating more processing of gas to stars. 
Our estimates of C,, computed assuming Mgas/Mdust = 
120, are presented in Table 5. We caution that the gas 
mass to dust mass ratio is highly uncertain. In SMGs, 
there is evidence suggesting it is « 60 (Kovacs et al] 
l2006f l. Adopting this lower value would imply lower 
gas masses by a factor of two and hence double our 
Mstar/Mgas estimates. 

Clustering studies suggest that the most luminous 
DOGs reside in ver y massive halos (Mdm ~ 10^^ Mq 
iBrodwin et al.lf2008l ) . It is tempting to attribute the low 
stellar masses we estimate for DOGs to youth. However, 
the absolute stellar masses we compute are extremely un- 
certain. For example, the use of a Salpeter IMF rather 
than a Chabrier IMF would approximately double our 
stellar mass estimates (Bruzual & Chariot 2003). Be- 
yond the choice of what IMF slope to use, the mass-to- 
light ratio of a model galaxy (for a given rest-frame near- 
UV - R color) can vary significantly depending on the 
details of its star-formation history, the dumpiness of its 
interstellar medium and the associated dust attenuation 
law, as well as how advanced stages of stellar evolution 
are treated, such as blue stragglers , thermally-pu lsating 
asymptotic giant branch stars, etc. (jConrov et al.l f20091. 
In light of these uncertainties, the fact that our stellar 
mass estimates are low (Mgtar 10^° Mq) compared to 
the dark matter haloes in which we believe they reside is 
not yet a cause for concern - a quantitative study of the 
maximum possible stellar mass allowed by the photom- 
etry (by examining results from different stellar popula- 
tion synthesis codes, star- format ion histories, metallici- 
ties, etc.) would be the best way to approach this issue 
in the near-term, but is beyond the scope of the current 
work. 

4. DISCUSSION 

In this section, we seek to understand the role of DOGs 
in galaxy evolution and their relation to other high-z 
galaxy populations. We begin by motivating a compari- 
son sample of such objects, including SMGs and Spitzer- 
selected ULIRGs from the XFLS and SWIRE survey We 
then examine how the measured properties differ from 
population to population. We end with the implications 
of these comparisons for models of galaxy evolution. 

4.1. Related z « 2 Galaxy Populations 



12 



Bussmann et al. 




Fig. 7. — contours based on SSP fits to HST imaging in H-band and V- or /-band. Black lines indicate 80%, 95%, and 99.9% 
confidence intervals. Grey contours trace lines of constant stellar mass for the given HST photometry, starting with A/star = 10^" Mq in 
the bottom left and increasing by 0.5 dex towards the upper right. In the top right corner of each panel is the minimum xi value and the 
associated stellar mass, in units of log(Mstar/A^Q). The bottom left corner contains the source identifier. Not shown are sources S2 and 
S7, since these targets have no HST imaging available. The best-fit stellar masses range from 10^0 - 10^2 Mq. 



4.1.1. SMGs 

SMGs represent an interesting population of galaxies 
for comparison with DOGs because they are selected 
at sub-mm wavelengths where the dominant emission 
component is cold dust (Tdust 30 K). In contrast, 
DOGs are selected predominantly by their brightness 
at 24/im and therefore should be dominated by hot 
dust. Despite this fundamental distinction, these two 
galaxy populations ha ve similar number densities and 
redshift distribut ions (IChapman et al.l I2005t IDev et all 
I2008t IBlain et al.ll2004f) . Recent evidence suggests that 
24^m-faint (^24^m ~ 0.1 — 0.5 mJy) DOGs have a 
composite SED whose shape in the far-IR closely mim- 
ics that of the av erage bright (Fg5o > 5 mJy) SMG 
(|Pope et al.ll2008af ). Futhermore, 24/.tm-faint DOGs and 
SMGs have similar real space correlation lengths (ro « 
6±2 Mpc), yet there is tentative evidence that DOG 
clustering strength increases with 24/Ltm flux density 



(rn » 13±3 h'^ Mp c for DOGs with F24^m > 0.6 mJy; 
IBrodwin et al.|[2008D . While these results are suggestive 
of an association between the two populations, the de- 
tails of such a connection are not yet clear. In an ef- 
fort to study this connection via their far-IR properties, 
we will compare the data presented in this paper with 
SHARC-II 350^m and MAMBO 1.2mm imaging of 25 
SMGs from the Submillimetrc Common User Bolome- 
ter Array (SC UBA) HAlf Degree Extragalactic Survey 
(iLaurent et al.l 2006: Kovacs et al.l I2006t iCoppin et al] 
l2008HGreve'et al..,2004l . 

4.1.2. XFLS Sources 

A set of Spitzer-selected galaxies from the 4 
deg^ XFLS sh are many proper ties with the 24^m- 
bright DOGs ()Yan et all [2001 . The specific se- 
lection criteria are similar, although not necessar- 
ily as extreme in their IR-optical flux density ratios: 
^24/.m > 0.9 mJy, j/Fj,(24/L(m)/j/F,,(8/Ltm) > 3.16, and 
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TABLE 5 

Dust Masses and Stellar Properties 







Afstar 




Source Name 


(10» Mq) 


(IQio Mq) 


Mstar/Mgas 


SI 


< 6.3 


(0.3 - 1.5) 


0.05 - 0.24 


S2 


7.9 ± 1.4 






S3 


4.9 ± 1.5 


< 10 


< 1.7 


S4 


< 5.5 


(1-3) 


> 0.15 - 0.45 


S5 


< 6.2 


(1.5-3) 


> 0.20 - 0.40 


S6 


< 7.5 


(10 - 20) 


> 1.1 - 2.2 


S7 


2.7 ±0.7 






S8 


2.5 ±0.9 


> 1 


> 0.33 


S9 


< 8.5 


(0.6 - 1.1) 


> 0.06 - 0.11 


SIO 


< 6.5 


(0.5 - 3) 


> 0.06 - 0.38 


Sll 


< 4.0 


(> 10) 


> 2.1 


S12 


< 1.9 


(1-2) 


> 0.44 - 0.88 



^ Dust mass assuming Tdust = 45 K. 
^ Stellar mass estimated from fitting photometry in 
the rest-frame UV and optical {V/I and H respec- 
tively). Range given reflects 95% confidence intervals 
based only on photometric uncertainty. 



i'F^{2Anm) /lyF^ {0.7 fim) > 10 (in comparison, DOGs 
have z/j;(24^m)/i/F^(0.7//m) > 30). Spitzer/IRS spec- 
troscopy of these objects has revealed strong silicate ab- 
sorption and in some cases PAH emission features on par 
with those of SMGs (fSaiina et al.ll2007f ). This suggests 
that the XFLS sources are composite AGN/starburst 
systems and may represent a trans ition phase betwee n 
(un)obscured quasars and SMGs ^ ma eTall \2{M ). 
MAMBO 1.2mm observations of 44 XFLS sources have 
allowed a detailed study of their fa r-IR prop erties and 
have suggested <iiR,> ~ TxlO^^L© ()Saiina et al...2008, ). 

4.1.3. SWIRE Sources 

The last set of comparison galaxies we consider 
are ^jpzteer-selected s ources from the SWIRE survey 
(jLonsdale et aLir2009( ). Like DOGs in Bootes and the 
XFLS sources, they have large IR-optical flux density 
ratios. However, an additional criterion has been ap- 
plied to identify sources with significant emission at rest- 
frame l.Gfim due to evolved stellar populations. For 
sources at z = 1.5 — 3, this means selecting objects whose 
mid-IR spectrum peaks at 5.8/.tm. Although spectro- 
scopic redshifts are not available for most of this sample, 
SED fitting has suggested photometric redshifts consis- 
tent with z ^ 2 and s tellar masses of (0.2 - 6)xl0^^ Mq 
(jLonsdale et al.|[2009t ). MAMBO 1.2 mm photometry for 
61 of these SWIRE sourc es has indicated far- IR lumi- 
nosities of 10i2_ioi3-3 (jLonsdale et al.ll2009l) . 

4.2. Comparison of Measured Properties 

Our results from section [3] represent our best estimates 
of LiR, ipiR, Tdust, Mdust, and M^tar for the DOGs in 
the sample. In TablelHl we give the median value of these 
quantities for DOGs in Bootes (from this paper), SMGs, 
and XFLS and SWIRE sources. In computing these me- 
dian values, we do not consider sources at z < 1; nor 
do we consider sources without detections at (sub-)mm 
wavelengths (see discussion on caveats to the analysis at 
the end of this section). Table [5] also makes a distinction 
between XFLS sources whose mid-IR spectra are dom- 
inated by strong PAH features (XFLS PAH) and those 



that show weak or absent PAH features (XFLS weak- 
PAH). Each of these galaxy populations is further sub- 
divided into those that quahfy as DOGs {R - [24] > 14) 
and those that do not. 

The primary feature of this comparison is that the rel- 
ative uncertainty in the estimated parameters between 
galaxy populations has been reduced by computing the 
respective values self-consistently with the methods out- 
lined in section [3l The exceptions to this rule are Lir 
and Tdust (note that while the photometry used to de- 
termine Mstar for Bootes DOGs is different than for the 
other galaxy populations, the methodology used is the 
same, including the use of the same set of model SSP 
templates). Our method of computing Ljr relies on the 
assumption that Mrk 231 represents a reasonable ap- 
proximation of the source SED. For many SMGs as well 
as XFLS and SWIRE sources, this is an unrealistic as- 
sumption. Instead, we estimate Lir from LpjR, assuming 
(1) T^iR = LiR_sB + TiR^ agn; (2) iiR, SB = aLpiR; (3) 
iiR, AGN = eLiR , where a is a fac tor w 1.3, depend- 
ing on Td and (3 (jHelou et al .1119881 ) and e is the typi- 
cal AGN fraction of the galaxy population. For SMGs, 
we adopt the conservative upper limit from iPope et al.l 
(jl008b) of e = 0.3. For SWIRE sources and XFLS 
PAH sources, i :his fraction is » 0.3 llPope et all l2008bl : 
ILonsdale et an[2009l : ISaiina et al.ll2008D. while for XFLS 
weak-PAH sources we use 0.7 (jSaiina etal|[2008l) . 

The Tdust values given in the literature are adopted 
for each source. It should be noted that Tdust for 
SWIRE sources are uncertain due to the lack of data 
near the far-IR peak (i.e., observed-frame 160 or 350 
^m). ILonsdale et al.l (|2009f ) analyze the stacked signal 
at 160/J,m from these sources and find that the Tdust is 
higher by as much as 10 K than what is assumed in their 
Table [6l For a given 1.2mm flux, increasing Td by 10 K 
will increase Tfir by a factor of «3 and decrease Mdust 
by «50%. 

The key result from Table [S] is that while the DOGs in 
our sample have lower dust masses than the other galaxy 
populations by a factor of '--^ 3 — 20, they have higher 
total IR and far-IR luminosities by factors of ~2. This 
distinction is driven by the difference in Tdust, as DOGs 
in the sample have higher values by «10-20 K compared 
to the other galaxy populations. 

In terms of the stellar and gas mass estimates, the 
relationship between DOGs in Bootes (i.e., the sample 
studied in this paper) and the remaining galaxy popula- 
tions is unclear. Even if a single dust temperature and 
a single dust to gas mass ratio for each of the sources 
studied in this paper is adopted, the uncertainties on the 
stellar mass estimates are large enough to allow greatly 
varying stellar mass to gas mass ratios. Sources satisfy- 
ing R — [24] > 14 (i.e., DOGs) tend to have higher gas 
masses compared to i? — [24] < 14 sources (under the 
assumption of a constant dust-to-gas mass ratio). This 
difference is at least in part due to a difference in dust 
temperatures; within this sample, DOGs have lower dust 
temperatures than non-DOGs. This is in contrast with 
the evidence for high dust temperatures in the Bootes 
DOGs studied in this paper and may be an indication 
that mm-detected DOGs represent a special subset of 
DOGs that is more representative of the mm-selected 
galaxy population than the DOG population. 

An important caveat to this comparison is that we are 
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TABLE 6 
Average High- 2 Galaxy Properties 









^IR 


^FIR 


'^dust 








Source 


i?-[24] 


TV 


(IOI2 Lq) 


(1012 Lq) 


(K) 


(109 Mq) 


(101° Mo) 




Bootes 


> 14 


5 


23 


10 


45 


0.5 


> 2 


> 0.3 


XFLS<^ 


ALL 


11 


7.7 


2.7 


32 


5.0 


13 


0.48 




^ 14 


p. 
u 


0. 1 




97 


( .0 


1 n 


21 






5 


6.5 


4.6 


37 


2.3 


16 


0.58 


XFLS PAH<: 


ALL 


5 


5.7 


1.8 


31 


4.6 


22 


0.54 




^> 


2 


yj.O 




9Q 


7 

/ .0 




n "^n 

U. JU 






Q 




9 n 


'i9 
oz 


9 s 
z.o 


22 


u.uy 


XFLS weak-PAH<= 


ALL 


6 


9.4 


3.4 


32 


5.3 


5.0 


0.25 




> 14 


4 


3.5 


0.8 


26 


7.2 


3.8 


0.17 




< 14 


2 


14 


8.0 


45 


1.6 


7.3 


0.40 


SWIRE^ 


ALL 


19 


6.5 


3.1 


32 


6.7 


28 


0.53 




> 14 


16 


6.7 


3.1 


32 


6.1 


28 


0.56 




< 14 


3 


5.5 


2.6 


30 


9.9 


24 


0.35 


SMG" 


ALL 


18 


6.9 


3.2 


35 


1.5 


10 


0.60 




> 14 


4 


7.8 


3.6 


28 


2.9 


6.9 


0.28 




< 14 


14 


6.1 


2.9 


37 


1.1 


12 


0.71 



^ Computed using Mgas/Mdust = 120. 

^ Includes only DOGs detected at 350/im. 

XFLS sources with MAMBO 1.2mm detections (|Saiina et alll2008l). 
SWIRE sources with MAMBO 1.2mm d etections (|Lonsdale et al.ll2009D . 

° From compilation of ICoppin et al.l (120081 ). 



dealing with small sample sizes due to incomplete cover- 
age at one or more bands from the mid- to the far-IR. 
For instance, while every DOG has a measured 24/im flux 
density, very few have been observed at 350/im, and only 
two have been observed at 1mm. Similarly, few XFLS 
and SWIRE sources have been detected at 1mm and even 
fewer have been observed at 350/im. Although SMGs are 
the best-studied class of objects within this set of popula- 
tions, they too suffer from low-number statistics. Larger 
sample sizes in the critical 200-500//m regime will arrive 
following the analysis of wide-field survey data from the 
Ballo on-borne Large Ap erture Submillimeter Telescope 
(e.g-. IPascale et al.ll2009 f) and the Herschel Space Obser- 
vatory. 

4.3. Implications for Models of Galaxy Evolution 

One of the major open questions in galaxy evolution 
is the effect that AGN have on their host galaxies. In 
the local universe, there is observational evidence that 
ULIRGs dominated by warm dust serve as a transition 
phase between cold dust ULIRGs and optically luminous 
quasars and that this transition may be driven by a major 
merger (|Sanders et al.ll988al lbl). Recent theoretical mod- 
els of quasar evolution based on numerical simulations of 
major mergers between gas-rich spirals have suggested 
that the subsequent growth of a super-massive black hole 
can regulate star-formation via a feedback effect which 
re-injects energy into the interstellar medium (ISM) and 
expels the remaining gas that would otherwise form stars 
(jHopkins et al.ll2006h . 

Although the notion that local UL IRGs are associated 
with mergers is well accepted (e.g., iSanders fc Mirabel 
Il996f ). the picture is less clear at high redshift. Mor- 
phological studies of high-z galaxies suffer from sur- 



face brightness dimming, making the detection o f faint 
merg e r remnant sign atures difficult (e.g., Dasvr a et al.l 
2008t iMelbourne et all 120081 : iBussmann et all 120091 : 
Melbourne et al.ll2009f l. However recent theoretical work 
on the cosmological role of mergers in the formation of 
quasars and spheroid galaxies suggests that they dom- 
inate the z > 1 quasar luminosity density compared 
to secular processes such as bars or disk instabilities 
(jHopkins et al.l[200l ). 

If major mergers drive the formation of massive galax- 
ies at high redshift, then one possible interpretation of 
our results involves an evolutionary scenario in which 
these sources represent a very brief but luminous episode 
of extreme AGN growth just prior to the quenching of 
star formation. In such a scenario, SMGs and the bright- 
est 24/im-selected sources represent the beginning and 
end stages, respectively, of the high star-formation rate, 
high IR luminosity phase in massive galaxy evolution. 
Consistent with this scenario is that we find 24/im-bright 
DOGs in Bootes to have higher dust temperatures (pos- 
sibly from AGN heating of the dust) than SMGs and less 
extreme Spiteer-selected sources. We caution, however, 
that these results are consistent with any evolutionary 
model in which the 24/im-bright phase follows the sub- 
mm bright phase, be it driven by major mergers, minor 
mergers, or some secular process. 

Finally, we stress that larger samples of mm and sub- 
mm imaging of Spzteer-selected galaxies are needed in 
order to understand their role in galaxy evolution fully 
by comparing samples of similar number density, clus- 
tering properties, etc. Much of this will be provided by 
upcoming Herschel and SCUBA-2 surveys. In the more 
immediate future, 1mm imaging with currently available 
instruments such as the Astronomical Thermal Emis- 
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sion Camera (AzTEC) and the MAx-Planck Millimeter 
Bolometer Array (MAMBO) will be critical to constrain- 
ing the cold dust properties of ^piteer-selected galaxies. 
Only when these surveys have obtained statistically sig- 
nificant numbers of detections or stringent upper limits 
will we be able to make definitive conclusions regarding 
the nature of the link between AGN and starbursts in 
the formation of the most massive galaxies. 

5. CONCLUSIONS 

We present CSO/SHARC-II 350^m and CARMA 1mm 
photometry of DOGs in the Bootes Field. The major 
results and conclusions from this study are the following: 

1. At 350/xm, 4/5 DOGs are detected in data with 
low rms levels (< 15 mJy) and 0/8 DOGs are de- 
tected in data with medium to high rms levels 
(20 — 50 mJy). At 1mm, a subset of two DOGs 
were observed but not detected. 

2. Mrk 231 is confirmed as a valid template for the 
SEDs of the DOGs in this sample. This suggests 
the 24/xm bright (i^24pm ^ 1 mJy) population of 
DOGs is dominated by warm dust, possibly heated 
by an AGN. Cold dust templates such as Arp 220 
are inconsistent with the data in all twelve objects 
studied. 

3. Trends in the flux density ratios 350/xm/24/xm 
and 1200/Ltm/24/im with the i?-[24] color 
(-F24/^m/-Fo.7/jm) show that DOGs in this sam- 
ple have elevated 24/xm emission relative to SMGs, 
most likely due to an obscured AGN. 

4. The non-detections at 1mm imply T^ust greater 
than 35-60 K for two objects. 

5. If the dust properties of the two DOGs observed 
at 1mm apply generally to the 24/im bright DOGs, 

then we estimate dust masses for these sources of 
1.6 — 6.1 X 10^ M0. Lower Idust would imply higher 
dust masses and vice versa. 

6. In comparison to other z « 2 ULIRGs, DOGs have 
warmer dust temperatures that imply higher IR lu- 
minosities and lower dust masses. This may be an 
indication that AGN growth has heated the ambi- 
ent ISM in these sources. 

7. Our stellar mass estimates provide weak evidence 
indicating that the 24/im-bright DOGs may have 
converted more gas into stars than SMGs or other 
Spitzer-selected sources, consistent with them rep- 
resenting a subsequent phase of evolution. An im- 
portant caveat to this conclusion is that we have as- 
sumed DOGs and SMGs share the same gas mass 



to dust mass ratio. Testing this assumption will 
require new data and will be an important goal of 
future work. 
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